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Abstract 

In this article, we study the vertices D*D*P, D*DV and DDV with the 
light-cone QCD sum rules. The strong coupling constants <?d*d*p, 9D*DPi 
Jd*dVi fD*D*v, 9ddv and gD*D*v play an important role in understanding 
the final-state interactions in the hadronic B decays. They relate to the basic 
parameters g, A and j3 in the heavy quark effective Lagrangian respectively. 
Our numerical values of the g, (3 and A are much smaller than most of the 
existing estimations. If the predictions from the light-cone QCD sum rules are 
robust, the final-state interaction effects maybe overestimated in the hadronic 
B decays. 
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1 Introduction 

Final-state interactions play an important role in the hadronic B decays [HI2]- How- 
ever, it is very difficult to take them into account in a systematic way due to the 
nonperturbative nature of the multi-particle dynamics. In practical calculations, we 
can resort to phenomenological models to outcome the difficulty. The one-particle- 
exchange model is typical (for example, one can consult Ref.[2]), in this picture, the 
soft interactions of the intermediate states in two-body channels with one-particle 
exchange make the main contributions. The phenomenological Lagrangian has many 
input parameters, which describe the strong couplings among the charmed mesons 
in the hadronic B decays [2]. In the following, we write down the relevant phe- 
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nomeno logical Lagrangian [2], 
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where we take the convention 60123 = 1- 

The strong coupling constants (for example, gn*DP, gD*D*p, etc.) can be esti- 
mated with the heavy quark effective theory and chiral symmetry [3]. In the heavy 
quark limit, the strong coupling constants in the phenomeno logical Lagrangian can 
be related to the basic parameters g, A and j3 in the heavy quark effective Lagrangian 
(one can consult Ref. [3 J for the heavy quark effective Lagrangian and relevant pa- 
rameters, here we neglect them for simplicity), 



gD*DP 2 

9d * d * p - TM^m -Jp 9 ' 
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where gv = 5.8 from the vector meson dominance theory [3]. For existing estimations 
of the values of the g, A and /3, one can consult Refs. |5lElCZllEll2llinilIIllI2llI3llI3llI5]. 

In previous work [15], we study the strong coupling constants of the DDV and 
D*DV with the light-cone QCD sum rules, the numerical values of the gnDV an d 
Jd*dv are much smaller than the existing estimations based on the vector meson 
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dominance theory [5]. In this article, we study the strong coupling constants go* d* p, 
Id'DVi 9ddv with the light-cone QCD sum rules 0. Furthermore, we analyze the 
corresponding parameters g, A, (3 in the heavy quark effective Lagrangian [3J, and 
obtain the values of the strong coupling constants gD*DP, fo*D*v and gn*D*v- 

The light-cone QCD sum rules carry out operator product expansion near the 
light 0, instead of short distance, x ~ 0, while the nonperturbative matrix 

elements are parameterized by the light-cone distribution amplitudes (which are 
classified according to their twists) instead of the vacuum condensates [HI [T7] . The 
nonperturbative parameters in the light-cone distribution amplitudes are calculated 
with the conventional QCD sum rules and the values are universal [18] . 

The article is arranged as: in Section 2, we derive the strong coupling constants 
gz>*D*p, Jd*dv an d Qddv with the light-cone QCD sum rules; in Section 3, the 
numerical result and discussion; and in Section 4, conclusion. 

2 Strong coupling constants gD*D*p, Id*dv an d Qddv 
with light-cone QCD sum rules 

We study the strong coupling constants gD*D*p, Jd*dv and gDDV with the two-point 
correlation functions IT^(j>, q), U^(p,q) and H^(p, q), respectively, 
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where the currents Jl{x) and J^{x) interpolate the mesons D*°, D* + , D*, D°, D + 
and D s , respectively. The i denote the u, d and s quarks respectively. The external 
states 7r, K, p, K* , (j) have the four momentum p^ with p 2 = m 2 , m\, m 2 p , m\», 
m\, respectively. 

According to the basic assumption of current-hadron duality in the QCD sum 
rules [18], we can insert a complete series of intermediate states with the same 
quantum numbers as the current operators J^x) and J${x) into the correlation 
functions H^ v (j), q), H^{p,q) and U^(p,q) to obtain the hadronic representation. 
After isolating the ground state contributions from the pole terms of the mesons D* 

2 In this article, we present the results for the strong coupling constants fo* dv and gDDV which 
are originally obtained in Ref.|15j explicitly, and perform a comprehensive analysis of the strong 
coupling constants in Eq.(l). 
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and Dj, we get the following results, 

fD*fD*M D *M D *g D * D * P 
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fD*fDM D *M 2 n f D * DV . 
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= 4(p, g )e F# y^ + -, (9) 

U lJ (p,q) = —p T 2e-q + --- 

(m + m c )( mj + m c ) {M 2 Di - (q + p) 2 } \M 2 D . - q 2 j 

= H%{p,q)e-q + --- , (10) 
where the following definitions for the weak decay constants have been used, 

f M 2 

(0|Jj(0)|Df(p)) = f D; M D .e„. (11) 

In Eqs.(8-10), we have not shown the contributions from the high resonances and 
continuum states explicitly as they are suppressed due to the double Borel trans- 
formation. Non-conservation of the vector currents J^x) can lead to non- vanishing 
couplings with the scalar mesons Dq, D + and D s0 , 

(0|Jj(0)|Ao(p)) = f DaPlt , (12) 

where the fu i0 are the weak decay constants. In this article, we choose the tensor 
structure € jJilia pp a q^ (or e liua pe l 'p a q 13 ) for analysis in Eqs.(8-9), the presence of the 
scalar mesons cannot result in contaminations. We have alternative choice to use 
the axial- vector currents ^(0) to interpolate the pseudoscalar mesons D°, D + and 
D s . However, the axial- vector currents J^(0) can also interpolate the axial- vector 
mesons D®, Df and D s i, 

(0|4(0)|Ai(p)> = f Dil M Dil e, : (13) 

where the fo u are the weak decay constants, we should be careful to avoid contam- 
inations from the axial-vector mesons. 

In the following, we briefly outline operator product expansion for the correla- 
tion functions n^(p, q), WHp,q) and IP J '(p, q) in perturbative QCD theory. The 
calculations are performed at large spacelike momentum regions {q + p) 2 <C and 
q 2 <C 0, which correspond to small light-cone distance x 2 ~ required by validity of 
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the operator product expansion. We write down the propagator of a massive quark 
in the external gluon field in the Fock-Schwinger gauge firstly |19j . 



The contributions proportional to G^ u can give rise to three-particle (and four- 
particle) meson distribution amplitudes with a gluon (or quark-antiquark pair) in 
addition to the two valence quarks, their corrections are usually not expected to 
play any significant roles. For examples, in the decay B — ► XcoK, the factorizable 
contribution is zero and the non-factorizable contributions from the soft hadronic 
matrix elements are too small to accommodate the experimental data [20]; the net 
contributions from the three-valence particle light-cone distribution amplitudes to 
the strong coupling constant Qd s1 d*k are rather small, about 20% [21]. In this 
article, we observe that the contributions from the three-particle (quark-antiquark- 
gluon) light-cone distribution amplitudes are less than 5% for the strong coupling 
constants go*D*p- The contributions of the three-particle (quark-antiquark-gluon) 
distribution amplitudes of the mesons are always of minor importance comparing 
with the two-particle (quark-antiquark) distribution amplitudes in the light-cone 
QCD sum rules. In our previous work, we also study the four form-factors fi{Q 2 ), 
f2(Q 2 ), 9\{Q 2 ) and (^(Q 2 ) of the S — ► n in the framework of the light-cone QCD 
sum rules approach up to twist-6 three-quark light-cone distribution amplitudes and 
obtain satisfactory results [22J. 

In a word, we can neglect the contributions from the valence gluons and make 
relatively rough estimations in the light-cone QCD sum rules. In this article, we 
take into account the three-particle light-cone distribution amplitudes of the pseu- 
doscalar mesons, and neglect the three-particle light-cone distribution amplitudes 
of the vector mesons to avoid cumbersome calculations. 

Substituting the above c quark propagator and the corresponding n, K, p, 
K* and mesons light-cone distribution amplitudes into the correlation functions 
H l i v (p, q), n*f (p, q) and ir j '(p, q), respectively, and completing the integrals over the 
variables x and k, finally we obtain the analytical results at the level of quark-gluon 
degrees of freedom. The explicit expressions are presented in appendix A. 

In calculation, the two-particle it, K, p, K*, <fi mesons and three-particle n, 
K mesons light-cone distribution amplitudes have been used JT6] [19], [23] [21], the 
explicit expressions are given in appendixes B-C. The parameters in the light-cone 
distribution amplitudes are scale dependent and can be estimated with the QCD 
sum rules [16] [19] [23] [21]. In this article, the energy scale p is chosen to be p c = 




(14) 
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Now we perform the double Borel transformation with respect to the variables 
Qi = —(.P + q) 2 an d Q\ = —q 2 for the correlation functions IT^., II?- and II?- in 
Eqs.(8-10), and obtain the analytical expressions of the invariant functions in the 
hadronic representation, 
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where we have not shown the contributions from the high resonances and continuum 
states explicitly for simplicity. 

In order to match the duality regions below the thresholds s® and s® for the 
interpolating currents, we can express the correlation functions Iljj, II? and Ilf- at 
the level of quark-gluon degrees of freedom into the following form, 
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{s t - (q + p) 2 } {s 2 - q 2 } 



(16) 



where the p^(si,s 2 



1,2,3) are spectral densities, then perform the double 
Borel transformation with respect to the variables Q 2 and Q\ directly. However, 
the analytical expressions of the spectral densities p%(si, S2) are hard to obtain, we 
have to resort to some approximations. As the contributions from the higher twist 

or -rU), the net contributions of 
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-{q+up) 2 



terms are suppressed by more powers of 
the twist-3 and twist-4 terms are of minor importance, less than 10% for the strong 
coupling constants Qd*d*Pi the continuum subtractions will not affect the results 
remarkably (for the strong coupling constants Gs(D s0 D*<p) and Ga(DsiDs4>)i the 
contributions are less than 20% [2S]-)- The dominating contributions come from 
the two-particle twist- 2 terms involving the 4>(u), <f>±(u) and <j>\\(u). We perform the 
same trick as Refs. [T9| [26] and expand the amplitudes 4>(u), <fi±(u), <p\\{u) and 4> a (u) 
in terms of polynomials of 1 — u, for example, 
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where the bk are coefficients, then introduce the variable s\ and the spectral densities 
are obtained. 

After straightforward calculations, we obtain the final expressions of the double 
Borel transformed correlation functions Iljj, II?- and 11^ at the level of quark-gluon 
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degrees of freedom. The masses of the charmed mesons are = 1.87GeV, Mo s 
1.97GeV, M D * = 2.010GeV and M D * = 2.112GeV, 
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0.87, 

0.78, (18) 
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there exist overlapping working windows for the two Borel parameters Mf and Mf, 
it is convenient to take the value Mf = Mf , 

Mf ~ Ml. ~ Ml. ~ Ml. ■ 1 J > 

We introduce the threshold parameters s — max^sf]) (s? and s° corresponding 
to M\ and Mf , respectively) and make the simple replacement, 

(\ 

m +un(l—un)m p m +uq(1 — iiQ)m p Sp 
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- e"^ (20) 

for the correlation functions Hjp n|- and Ilf- respectively to subtract the contribu- 
tions from the high resonances and continuum states [19] . 

Finally we obtain the sum rules for the strong coupling constants Qd*d*Pi ]d*dv 
and g DD v, 
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Here we write down only the analytical results without the technical details, one 
can consult appendix D for some technical details. 

In the following, we present another approach for subtracting the contributions 
from the high resonances and continuum states [2H EI]- Firstly, we perform a 
double Borel transformation with respect to the variables Q\ and Q\ respectively, 
and obtain the result, 
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where f{u) stand for the light-cone distribution amplitudes, Aj = A 2 = m 2 + M (l — 
w )p 2 and p(si, s 2 ) stand for the corresponding spectral densities. Then, we make a 
replacement Mf — > — , M| — > — in above equation, 
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Finally, we take a double Borel transformation with respect to the variables o\ and 
cr 2 respectively, the resulting QCD spectral density reads 



i.e. 
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For the twist-2 terms, a = 1, the two subtracting approaches lead to the same 
results, the simple subtraction procedure we take in Eq.(20) is still reasonable. 



3 Numerical result and discussion 

The input parameters are taken as m s = (0.14 ± 0.01)GeV, m c = (1.35 ± 0.10)GeV, 
m u = md = (0.0056 ± 0.0016)GeV, f K = 0.160GeV, f w = 0.130GeV, f p = 
(0.216 ± 0.003)GeV, = (0.165 ± 0.009)GeV, f K , = (0.220 ± 0.005)GeV, = 
(0.185 ±0.010)GeV, fa = (0.215 ± 0.005)GeV, f£ = (0.186 ±0.009)GeV [24J, m K = 
0.498GeV, = 0.138GeV, m p = 0.775GeV, m K * = 0.892GeV, m = 1.02GeV , 
M D = 1.87GeV, M Ds = 1.97GeV, M D * = 2.010GeV and M D * = 2.112GeV. 

For the K meson: A 3 = 1.6±0.4, f 3K = (0.45±0.15) x l(T 2 GeV 2 , lu 3 = -1.2±0.7, 
u 4 = 0.2 ± 0.1, a 2 = 0.25 ± 0.15, a x = 0.06 ± 0.03, 774 = 0.6 ± 0.2 pH HI [23]. 

For the ir meson: A 3 = 0.0, / 37r = (0.45 ± 0.15) x 10~ 2 GeV 2 , uj 3 = -1.5 ± 0.7, 
u 4 = 0.2 ± 0.1, a 2 = 0.25 ± 0.15, a x = 0.0, rj A = 10.0 ± 3.0 [TBI [HI [23]. 

For the p meson: a\ = 0.0, af = 0.0, a| = 0.15 ± 0.07, = 0.14 ± 0.06, 
Cj = 0.030 ± 0.010, A3 1 = 0.0, £ij = -0.09 ± 0.03, 4 = 0.0, dj = 0.15 ± 0.05, 
A} = 0.0, 4 = 0.0, uj^ = 0.55 ± 0.25, = 0.0, C4 = 0.15 ± 0.10, £[ = 0.10 ± 0.05 
and (I = -0.10 ±0.05 [21]. 

For the K* meson: of = 0.03 ± 0.02, a{ = 0.04 ± 0.03, a\ = 0.11 ± 0.09, 

4 = 0.10 ± 0.08, (I = 0.023 ± 0.008, aJ = 0.035 ± 0.015, = -0.07 ± 0.03, 
4 = 0.000 ± 0.001, ujI = 0.10 ± 0.04, A3 1 = -0.008 ± 0.004, nj = 0.003 ± 0.003, 
i4 = 0.3 ± 0.1, A^ = -0.025 ± 0.020, C4 = 0.15 ± 0.10, Cj = 0.10 ± 0.05 and 
CJ = -0.10 ±0.05 [21]. 

For the (p meson: a\ = 0.0, af = 0.0, o| = 0.18 ± 0.08, = 0.14 ± 0.07, 
(J = 0.024 ± 0.008, aJ = 0.0, uj\ = -0.045 ± 0.015, 4 = 0.0, Jj = 0.09 ± 0.03, 
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M gr (GeV) 


M gr (GeV)(exp) 


v^(GeV) 


M 2S (GeV) 


O-(cn) 


1.90 ± 0.03 


1.869 


2.45 ± 0.15 


2.589 


l~(cn) 


2.00 ± 0.02 


2.010 


2.55 ± 0.05 


2.692 


O-(cs) 


1.94 ± 0.03 


1.969 


2.50 ± 0.20 


2.700 


l-(cS) 


2.05 ± 0.04 


2.112 


2.65 ± 0.15 


2.806 



Table 1: Numerical values of the ground state masses M gr and the threshold pa- 
rameters y/so from the QCD sum rules [32]. We denote the first excited state as 2S 
state, the values of the masses of the 2S states are taken from the predictions of the 
quark model [33] . 





9d*d*p 


9ddv 


!d*dv 


aS(GeV) 


6.5 ±0.5 






^(GeV) 


7.0 ±0.5 






s»(GeV<) 




6.0 ±0.5 


6.5 ±0.5 






6.3 ±0.5 


7.0 ±0.5 


s»(GeV<) 




6.3 ±0.5 


7.0 ±0.5 



Table 2: Threshold parameters for the strong coupling constants (?d*d*p, 9ddv and 
Id*dv, respectively. 

Aj = 0.0, 4 = 0.0, u£ = 0.20 ± 0.08, = 0.0, C4 = 0.15 ± 0.10, (I = 0.10 ± 0.05 



and d = -0.10 ± 0.05 

The values of the decay constants fa, fr> s , fp* and fo* vary in a large range 
from different approaches, for example, the potential model, QCD sum rules and 
lattice QCD, etc [271 [28]. F° r the decay constant fp, we take the experimental data 
from the CLEO Collaboration, f D = (0.223 ± 0.017)GeV [29]. If we take the value 
f Ds = (0.274±0.013)GeV from the CLEO Collaboration, the 577(3) breaking effect is 
rather large, 4^ = 1.23, while most of the theoretical estimations indicate 4^ ~ 1.1. 

JD JD 

In this article, we take the value = 1.1. For the decay constants fo* and fr>*, 

JD ' ' " s 

we take the central values from lattice simulation [3D], /d* = (0.23 ± 0.02)GeV and 
f D * = (0.25 ± 0.02)GeV, 

LEI ~ ^ = 1.1. (29) 

JD* JD 

The duality threshold parameters So are shown in Table. 2, the numerical (central) 
values of So are taken from the QCD sum rules for the masses of the pseudoscalar 
mesons D°, D + , D s and vector mesons D*°, D* + , D*, see Table. 1 [32J. The threshold 
parameters and s° (corresponding to M\ and M| respectively) are not equal 
in some channels in Eqs. (21-23), we choose the larger one i.e. sq = max^^s!]) 
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to take into account all the contributions from the ground states, certainly, there 
maybe some contaminations from the 2S state in the channel with smaller threshold 
parameter i.e. min^, s^), and impair the predictive power. The uncertainties of 
the threshold parameters sq are about Ssq = (0.25 — 1.0)GeV 2 , see Table. 1, in this 
article, we take 5s = 0.5GeV 2 for simplicity. 



The Borel parameters are chosen as M{ 



M 2 and M 2 



(3 - 7)GeV 2 , in those 



regions, the values of the strong coupling constants go*D*p, 9ddv and Jd*dv are 
rather stable. 

In the limit of large Borel parameter M 2 , the strong coupling constants gD*D*p, 
fo*DV and gnov take up the following behaviors, 
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M 2 f Vlj a\ 
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(30) 



It is not unexpected, the contributions from the twist- 2 light-cone distribution am- 
plitudes 4>(u), 4>\\(u) and (f>±(u) are greatly enhanced by the large Borel parameter 
M 2 , uncertainties of the relevant parameters presented in above equations have sig- 
nificant impact on the numerical results. 

Taking into account all the uncertainties, finally we obtain the numerical values 
of the strong coupling constants gD*D*p, Jd* dv and gDDV, which are shown in 
Figs. 1-3, respectively, 



gD*D*n 
gD*D*K 

Id* Dp 

fo*D s K* 
fD*D s <f> 
gDDp 

9dd s k* 

gD 3 D s( p 

The average values are about 

gD*D*p 
Id*dv 

gDDV 



(3.30 ± 1.55)GeV _1 
(3.50 ± 1.57)GeV _1 
(0.89±0.15)GeV- 1 
(1.01±0.20)GeV _1 
(0.82±0.16)GeV" 1 
1.31 ±0.29, 
1.61 ±0.32, 
1.45 ±0.34. 



(3.40 ± 1.55)GeV _1 
(0.91 ±0.17)GeV" 1 
1.46 ±0.32. 



(31) 



(32) 



The corresponding values of the basic parameters g, A and f3 in the heavy quark 
effective Lagrangian can be obtained from Eq.(3), and listed in Tables. 3-5. From 
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Tables. 3-5 or Eq.(3), we can obtain the values of the strong coupling constants 

9D*DPi fo*D*V and gD*D*Vi 

9d*dp = 6.73 ±3.07, 

f D *D*v = 1.85 ±0.35, 

g D * D *v = 1.46 ±0.32. (33) 

Taking the replacements goD P - > 2£2r an d fn*D P - > ^ D *°p i n Eq.(l), we can obtain 
the same definitions for the strong coupling constants in Ref. |31j. Our numerical 
values gDDp = 2.62 ± 0.58 and fo*D P = (3.56 ± 0.60)GeV _1 are compatible with the 
predictions g DDp = 3.81±0.88 and f D * Dp = (4.17±1.04)GeV _1 in Ref. [31]. The value 
gD*DP = 6.73 ± 3.07 obtained from the relation of heavy quark effective theory is 
different from the one obtained with the light-cone QCD sum rules gn* on = 12.5±1.0 
[j~9] . In Ref. [31] . the authors take much smaller values for the decay constants of 
the charmed mesons than the present work. In Ref. [19], the authors take much 
smaller value of the decay constant fo and much larger value of the nonperturbative 
parameter 02 (/^) than the present work. It is not unexpected that the numerical 
values are different from each other, see Eq.(30). We can expect the relations in 
Eq.(3) work well. 

The values of the g vary in a large range from different approaches, see Table. 3, 
the present prediction g = 0.22 ± 0.10 is consistent with our previous calculation 
g = 0.16^0 05 with the light-cone QCD sum rules [II]. However, it is much smaller 
than most of the existing estimations, this maybe due to the shortcomings of the 
light-cone QCD sum rules. 

The basic parameter A relates to the form-factor V (q 2 ) of the hadronic transitions 
(V I <Z7a»(1 — 75)^ I B) and (V \ qcr^il + 75)6 | B), which can be calculated with 
the light-cone sum rules approach and lattice QCD. With assumption of the form- 
factor V(q 2 ) at q 2 = q^ ax = (Mb — My) 2 is dominated by the nearest low-lying 
vector meson pole, we can obtain the values of the A H|, which are presented 
in Table. 4. The parameter (3 can be estimated with the vector meson dominance 
theory, which is presented in Table. 5, for technical details, one consult Ref. [15] . The 
large discrepancies maybe that the vector meson dominance theory overestimates 
the values of the [3gy and Xgy, the other possibility maybe the shortcomings of the 
light-cone QCD sum rules. 

We can borrow some idea from the strong coupling constant go*D-Ki the central 
value (go'D-w = 12.5 or go*D-K = 10.5 with the radiative corrections are included 
in) from the light-cone QCD sum rules is too small to take into account the value 
(gD*Dn = 17.9) from the experimental data [Hi [19j 154] . Naively, we can expect 
that the contributions from the radiative corrections cannot smear the discrepancies 
between our predictions and other estimations for the strong coupling constants 
gD*D*p, Jd*dv and gDDV- It has been noted that the simple quark-hadron duality 
ansatz which works in the one-variable dispersion relation might be too crude for 
the double dispersion relation [35]. As in Ref. [31], we can postpone the threshold 
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M 


Reference 


0.38 ±0.08 


[3J 


0.34 ±0.10 


[7J 


0.28 


[8J 


0.35 ±0.10 


[9J 


0.50 ±0.02 


[10] 


0.6 ±0.1 


[HJ 


0.59 ±0.07 


[12j 


n 97+U.U6 

u - z ' -0.03 


[13J 


lfi+ u - UY 


m 


0.22 ±0.10 


This work 



Table 3: Numerical values of the parameter g. 



(1 


Reference 


0.9 





0.36 ±0.08 


This work 



Table 4: Numerical values of the parameter (3. 



|A|(GeV-) 


Reference 


0.56 


[5J 


0.63 ±0.17 


[6J 


0.22 ±0.04 


This work 



Table 5: Numerical values of the parameter A. 
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— * — Central value; 






• Upper bound; 






A Lower bound. 
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M 2 (GeV 2 ) 



Figure 1: gD*D*n{A) and gD*D* a K{^>) with the Borel parameter M 2 after taking into 
account all the uncertainties. 

parameters Sq to larger values to include the contributions from the radial excitations 
(D' or D*') to the hadronic spectral densities, with additional assumption for the 
values of the gD*'D*Pi 9d'dv, fo*'DVi Id*d'Vi etc, we can improve the values of the 
gD*D*p, 9ddv an d fD*DVi an d smear the discrepancies between our values and other 
predictions. It is somewhat of fine-tuning. 

From Tables. 3-5, we can see that our numerical values are much smaller than 
most of the existing estimations (for example, the values taken in Ref.[2]). Naively, 
we can expect that smaller values of the strong coupling constants gD*D*p, 9d*dp, 
Id*dv, Id*d*v, Qddv and gD*D*v lead to smaller final-state interaction effects in the 
hadronic B decays. For example, the contributions from the rescattering mechanism 
in the decay 

B -> D*p -> Dtx 

can occur through exchange of D* (or D) in the t channel for the sub-precess D*p — ► 
Dn [2j. The amplitude of the rescattering Feynman diagrams is proportional to 

C 1 g D * D ^fD*Dp + C 2 g D *DngDD P oc D^X + D 2 gP , (34) 
where the Cj and are some coefficients. 



4 Conclusion 

In this article, we study the vertices D*D*P, D*DV and DDV with the light-cone 
QCD sum rules. The strong coupling constants g D *D*p, gD*DP, Id*dv, Id*d*v, gDDV 
and gn*D*v play an important role in understanding the final-state interactions in 
the hadronic B decays. They relate to the basic parameters g, A and (3 respectively 
in the heavy quark effective Lagrangian. Our numerical values of the g, (5 and A 
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Figure 2: g£>£> p (A), Qdd 3 k* (B) and gD s D a if>{C) with the Borel parameter M 2 after 
taking into account all the uncertainties. 
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Figure 3: /o*Dp(A), /d*d s /^*(B) and fD*D s ct>(C) with the Borel parameter M 2 after 
taking into account all the uncertainties. 
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are much smaller than most of the existing estimations. If the predictions from 
the light-cone QCD sum rules are robust, the final-state interaction effects maybe 
overestimated in the hadronic B decays. 



Appendix A 

The explicit expressions of the ITJ-, U 2 j and IT^- at the level of quark-gluon degrees 
of freedom, 



n 



n 



where 



u 



AA 3(mj + irij 



3) JO 



du- 



AA 2 



duA{u) 



1 2m 2 
AA 2 + AA 3 ~ 



~ fp ^ m k J Q da a J o da i J o dv J^ 



\u=ctj+(l— v)a g 



V\\ v± 



1 oij a g ,a g ,a>j) 



i r i 



1—a.a 



-fp,3 m P l3 / dv / da 9 

JO Jo •/() 

A 4 A 4 ^ 4- Vi 



da. 



Jo 



da 



d 1 



r 1 r 1 r a 9 /-1-/3 

+/p lJ m p li / / da 3 df3 da 
Jo Jo Jo Jo 



du AA 2 l^+a-*)"* 

1 — a — a 9 , a) 

d 1 



|m=1— va a 



[1 — a — (5, (3, a) + 



(35) 



+ 



AA 



± rrij + mj 
m Vi . 



o 



1 



+ 



2ml 



AA 2 AA 3 



'"'■"' v » I' 1 ^ 9± ) ( u ) 
mi + mj 



Jo 



(36) 



duA(it) 



m c m v .. I du 



hf{u) 
AA 2 



1 



+ 



2m 2 . 



7 rfr / 

Jo JO 



AA 2 AA 3 
dtC(t) 



1 2/// 2 



AA 2 AA 3 



+ 



(37) 



AA = m 2 c — (q + upY 
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Appendix B 

The light-cone distribution amplitudes of the K meson are defined by 

/l s 2 2 

2 p ■ x J 

(0\u(0)i^s(x)\K{p)) = jKW?K f due-^Mu; , 

m s + m u Jo 



(0\u(0)W5s(x)\K(p)) = i{p„x v -p v x») fKV f K f due- iup - x 4> a (u) , 

OyTTlg + fTl u ) 







(0\u(0)a^ 5 g s G af3 (vx)s(x)\K(p)) = f 3K {(p^p a g^ - PvP a g^) - (p^PpgL 

(0\u(0)-f ll -f 5 g s G a fs(vx)s(x)\K(p)} = f K m 2 K pJ^-^ — 

p • x 

J V ai A {l (a^e- ip - x(as+vas) 
+fKm 2 K (p^-p a g^) 



(0\u(0)-f lx ig s G a/3 (vx)s(x)\K(p)) = fKm 2 K p, ? aXp - V J Xa 

2 



+fKm K ( Pf3 g afi -p a g^) 
J Va t V ± (a t )e- tp - x{as+va3) , (38) 



where g^ = 9iXV - v ^ v x vX \ = ^e^G^ and Vol, = da u da s da 9 5(l 

^ - Oig). 



a, 
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The light-cone distribution amplitudes of the K meson are parameterized 



4>(u,ii) 
i> p (u,fi) 



<f>3(ati,fi) 



A ± (a i ,fj l ) 
A(u,fi) 



g(u,(i) 
B(u,fi) 



6u(l-u){l + ai cl(C) + a 2 cl(0} , 
1+{30^-^ 2 }c|(O 



+ 

6m(1 — u 



1 7 2 3 2 

360a u a s a 2 g jl + A 3 (a M - a s ) + u 3 ^(7a g - 3) 

120a u a s a g (v 00 + Vi (3a 9 - 1)) , 

120a u a s a g a 10 (a s - a u ) , 

-30a 2 {h 00 (l - a g ) + hoi [a g (l - a g ) - Ga u a s 



c|(0 



10 



30a 2 (a u - a s ) j/i o + V" 9 + ^h 10 (5a g - 3)| 



6m(1 — u 
+ 



16 24 20 

is + 35 a2 + 20m + y 774 



117 10 

"15 + 16 " 27^ 3 " 27* 



+ 



n 

210 



&2 



4 

135 



V3^3 



18 



-a 2 + 217740^4 



{2m 3 (10 - 15m + 6m 2 ) logM + 2m 3 (10 - 15m + 6m 2 ) log 
+mm(2 + 13mm)} , 

g(u,n) - <j)(u,fi) , 



u 
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where 

u - ^ 
n oo — v oo — — y , 

21 9 
«io = y ^4^4 - 7^a 2 , 

21 

o 

7 7 3 

ftoi = ^4^4 - — a 2 , 

7 3 
"io = t^ 4 ^ 4 + 20 a2 ' 

18 20 
#2 = 1 + y a 2 + 60% + y ?7 4 , 

9 

9a = -^2-6773^3, (40) 

here £ = 2« - 1, and Cf (£), CJ(£), (£), cf (£), cf (0 are Gegenbauer polyno- 
mials, 773 = ipuhi+iik an d p 2 = K+^) a [TBI [13 [23]. The corresponding light-cone 
distribution amplitudes for the n meson can be obtained with a simple replacement 
of the nonperturbative parameters. 
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Appendix C 



The light-cone distribution amplitudes of the K* meson are defined by 

-1 ( ^2 J2 



(0\u(0) llx s(x)\K*(p)) = Pli f K .m K . 



(0\u(0)s(x)\K*(p)) 
\u(0)a^s(x)\K*(p)) 



(0\u(0)w 5 s(x)\K*(p)) 



+ 



€ ■ X 

p-x J n 
e ■ x 



due~ mpx U^u) + T ^^A{u) 



16 



1 

-x 



p ■ x 
e ■ x 



f K ,m K , f due-^ x g K l\u) 
Jo 

* l l ..... 



'■xAv), 



v Jk^ k , / due-^ x C(u) 

* [p-x) J 



fx* ~ Ik 



m 2 K *e-x I due- lup - x hf(u) , 
'o 

2 2 

m K *x 



e ■ x 



(p-xj Jo 



f 

Jo 



16 

due~ iup - x B , (u) 



P-x Jo 



due- iup - x C ± (u) , 



Ik* - f K * 



m K * 



m K *e f t Ua pe u p a x 13 



due- lup - x g { l\u) 



(41) 
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The light-cone distribution amplitudes of the K* meson are parameterized as 



<j>\\{u,n) = 6 M (l- M )|l + a!3e + a^(5e 2 -l)| , 
<l>±(u,») = 6 M (l- M )(l + a^ + o^(5e 2 -l) 



<£W) = ^(l + a + a!^ 3 + {^5 + 5C|}(3e 2 -l) 
+ {54-^Al + ^A|}e(5e 2 -3) 

+ {^4 + ^I-^S}(3-30e 2 + 35a, 
<7?W) = 6u«|l+Qo™ + ^«!)c7 1 i (0+ 

G a2 + y cS + ^ - 2i" s ) ^(o + (i^S - 1>8) cm } 

/»<*W) = 3e 2 + ^(3e 2 -i) + ^4e 2 (5e 2 -3) + ^(3-3oe 2 + 35e 4 ) 

+ (y4 " ^) e(5? 2 - 3) , 

= i + {-i-^I + fcl-fc 4 }c|(0 

+ r28 a2+ 4 C3 "l6^"Y a;3 / Q(0 ' 

hs(u,fi) = i + {-i+^4-io(c 4 T +cr)}c|(o + {-^-^}ci(o, 

on „f4 4 || 8 ji 20 . 
= 30 W 2 m 2 |- + — 4 + + 

C(u,fi) = g 3 (u,n) + <j>\\(u,n) - 2g { l\u,/j,) , 

1 1 

B±(u,fi) = h§\u,n) - -(f) ± (u,ii) - -h 3 (u,[i) , 

C±(u,fj,) = h 3 (u, ji) - <f)j_(u, ji) , (42) 

where £ = 2u - 1, and CJ(£), C/(0, Cf (0, C 2 f (£), Cf(0 are Gegenbauer polyno- 
mials. The corresponding light-cone distribution amplitudes for the p and <f> mesons 
can be obtained with a simple replacement of the nonperturbative parameters. 
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Appendix D 



Here we present some technical details necessary in performing the Borel transfor- 
mation which are not familiar to the novices, 



d 

OLg) — exp 
du 



m c + u(l — u)m 



M 2 



r 1 r 1 r 1 -^ 

dv da g / da s f(v,a s , 
Jo Jo Jo 

S(U - U )\ u =a s +(l~v)a g 

rl rl rl fl-oig 

= I du I dv I da g / da s f(v, a s , a g )8 [u — a s — (1 — v)a g \ 
Jo Jo Jo Jo 

m 2 + u(l — u)m 2 K 



— exp 

du 



M 2 



8(u - U ) 



pi ru rl—a s 

= I du I da s / da. 

JO JO Ju-a s 



f(v,a s ,a g ) d 
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a g du 

2 i ,,/ \ „.\^2 

c 



exp 



m 2 + u(\ — u)m 2 K 



M 2 



m 2 + w(l — u)m 2 K 



M 2 



d r 



S(u - u o)-j- I da s 



l—a s 



da, 



S(u - u ) 
f(v,a s ,a g ) 



a n 



= — exp 



ml + u (l - u )m 2 K 
M 2 



d 

dU J:, 



da* 



l-a s 



da, 



f(v,a s ,a g ) 



where the f(v, a si a g ) stand for the three-particle light-cone distribution amplitudes. 
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